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REMOTE  SENSING  OF  ATMOSPHERE  VISIBILITY:  A  CRITICAL  REVIEW 


1.  INTRODUCTION 
1 .  1  Purpose 

This  is  a  final  report  on  Contract  F19628-77 -C -0202,  Remote 
Visibility  Study,  prepared  for  Air  Force  Geophysics  Laboratory  (OPA). 

The  purpose  of  this  study  was  to  conduct  a  critical  review  of  the  published 
concepts  for  satellite'  based  remote  probing  of  the  optical  properties  of 
the  lower  troposphere  particularly  the  aerosol  properties.  Further  the 
study  was  to  assess  the  applicability  of  these  and  any  other  methods  to  the 
determination  of  slant  path  contrast  reduction  and  atmospheric  transmission 
for  aircraft  viewing  geometries.  These  two  factors  control  the  ability 
of  optical  and  electro-optical  systems  to  detect  or  sense  targets. 

Transmission  affects  the  amount  of  light  which  reaches  the  system  from 
the  target  while  contrast  reduction  changes  the  ratio  of  light  from  the  target 
to  light  from  the  surrounding  background,  i.e.  the  rargel  contrast.  Whim 
contrast  falls  below  a  threshold  or  minimum  value,  the  system  can  no  longer 
detect  or  see  the  target.  The  major  contrast  reducing  agent  is  the  aerosol 
in  thi>  atmospheric  which  scatters  lighl  into  the  system. 

It  has  become  common-place  to  use  the  word  visibility  to  de¬ 
scribe  the  atmospheric  state  as  it  affects  the  target  detection  capability 
of  an  electro-optical  system  even  when  the  visible  part  of  the  optical  spectrum 
is  not  involved.  The  less  common  term  "see-ability"  is  more  precise  but 
still  is  misleading  by  implying  a  restriction  to  the  visible  spectrum.  With 
the  explicit  recognition  of  its  drawbacks,  the  word  visibility  will  be  used 
in  this  report  for  convenience. 

A  need  for  close  aerial  reconnaissance  exists  even  in  this  age 
of  globe -encircling  reconnaissance  satellites.  Since  the  reconnaissance 
plane  is  put  into  considerable  risk  in  executing  its  mission,  it  is  important 
to  have  prior  knowledge  that  the  technical  factors  are  favorable.  The  most 


important  factor  is  visibility  conditions  over  the  target.  The  military 
commander  would  have  an  important  tool  if  he  had  the  means  of  remotely 
sensing  visibility.  The  artificial  earth  satellite  is  an  obvious  candidate  for 
this  purpose. 

Indeed,  satellites  are  used  today  in  other  remote  sensing 
roles.  With  some  degree  of  success,  they  are  used  to  determine  the 
vertical  air  temperature  profile  for  meteorological  uses.  Other  uses 
include  remote  assessment  of  geological  and  agricultural  surface  conditions. 
Examination  of  the  imagery  from  the  latter  applications  (Ref.  1)  reveals 
that  the  presence  of  the  earth's  atmosphere  distorts  measurements  of  the 
radiance  of  surface  objects  and  corrections  must  be  applied  to  arrive  at  the 
true  radiance.  Since  these  atmospheric  effects  are  the  same  that  produce 
changes  in  reconnaissance  visibility,  there  is  some  basis  to  the  conjecture 
that  satellite  observations  can  be  used  to  deduce  surface  visibilities.  How¬ 
ever,  there  is  no  direct  path  from  these  qualitative  observations  to  the 
quantitative  methods  of  measurement  that  are  sought. 

The  concept  of  remote  visibility  sensing  from  satellites  has  been 
suggested  and  explored  earlier  by  numerous  investigators.  Their  published 
work  has  been  reviewed  and  the  most  promising  methods  selected  for 
presentation.  Some  new  concepts  which  arose  during  the  course  of  this  review 
are  also  presented.  Brief  mention  is  made  of  some  rejected  methods. 

1 .  2  Summary 

The  significant  and  practical  literature  on  remote  visibility 
sensing  from  satellites  reviewed  in  this  report  covers  the  interval  1970  to 
1979.  Of  course  no  experimental  literature  on  this  subject  could  have  existed 
prior  to  1958.  However,  the  idea  of  using  artificial  earth  satellites  to 
examine  the  earth's  atmosphere  was  beginning  to  be  considered  prior  to  Sputnik 
as  evidenced  by  Reference  2  which  proposes  a  method  of  determing  the 


vertical  ozone  distribution.  Still  by  1972,  little  work  on  visibility  sensing 
had  been  published.  An  extensive  review  article  on  the  general  problem 
of  atmospheric  sensing  from  satellites  covering  up  to  1972  (Reference  3) 
cited  only  four  archival  literature  references  related  to  particulate  sensing. 

The  literature  search  was  based  on  an  earlier  bibliography  prepared  for 
AFGL  (Reference  4)  and  was  extended  by  a  search  in  the  major  English 
language  optical,  meteorological,  and  environmental  journals.  Also  included 
in  this  review  was  a  series  of  international  conferences  on  remote  sensing 
in  the  environment. 

Remote  sensing  of  atmospheric  visibility  reduces  essentially 
to  detecting  the  amount  of  particulate  matter  in  the  first  5  kilometers  of 
the  troposphere.  Methods  can  be  characterized  in  one  way  by  whether  the 
illuminating  source  is  active  or  passive.  Active  methods  use  lidar  systems 
and  their  use  over  international  regions  might  bo  restricted  bv  political 
considerations.  The  applicability  of  active  methods  is  summarized  from 
more  extensive  studies  of  other  researchers.  The  most  important  passive 
light  source  is  the  sun  which  can  be  viewed  directly,  by  specular  reflection 
from  water  surfaces,  or  by  diffuse  re  flection  from  atmosphere  and  earth's 
surface  together.  Probing  the  atmosphere  by  viewing  the  sun  directly 
through  the  atmospheric  limb  of  the  earth  has  only  succeeded  in  reaching 
levels  about  10  km  above  the  surfaces  (References  5  and  6).  Thus,  this 
concept  is  not  applicable  to  this  problem.  Viewing  the  sun  through  specular 
reflection  off  of  a  water  surface  has  been  proposed  as  a  method  of  measuring 
atmospheric  transmission  but  no  tests  of  the  technique  have  been  made.  Most 
of  the  methods  which  have  been  tested  with  actual  satellite  data  arc'  based 
on  measurements  of  the  upwelling  light  diffusely  reflected  from  the  earth.  This 
method  is  the  easiest  to  implement  but  is  limited  to  use  over  surfaces  of  low 
reflectance ,  in  particular  the  ocean.  Any  method  which  works  for  the  sun  could 
be  used  with  the  moon  as  source  as  long  as  the-  measuring  instruments  have 


sufficient  sensitivity.  Another  possible  passive  source  is  the  night 
illumination  from  cities.  A  method  using  this  source  is  developed  from 
a  modification  of  the  Langley  method  of  measuring  the  atmospheric  trans¬ 
mission  using  rhe  sun. 

1'he  report  is  organized  as  follows:  first,  a  brief  discussion 
of  the  technical  background  sufficient  to  emphasize  the  major  factors  and 
to  define  various  terms  is  presented.  Then  follows  the  technical  discussions 
of  the  various  promising  techniques.  Some  rejected  concepts  are  briefly 
mentioned.  Finally,  the  conclusions  and  recommendations  derived  from 
this  study  are  presented. 
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2.  1  Atmospheric  Properties 

The  visual  appearance  of  the  atmosphere  is  so  familiar  and 
such  a  common  experience  that  we  seldom  inquire  deeply  as  to  '.he  causes 
of  the  ever  changing  panorama  presented  to  our  eyes.  Thi  '  '  to  see 

distant  objects  clearly  (sometimes  imprecisely  called  visibility)  changes 
from  excellent  on  those  brilliant  clear  blue  sky  days  following  the  passage 
of  a  cold  front  to  poor  on  those  hot,  humid,  and  hazy  summer  days.  These 
changes  in  visibility  are  controlled  by  changes  in  atmospheric  properties 
which  will  be  described  in  this  section.  This  material  is  given  to  provide 
the  reader  with  the  minimum  background  needed  to  read  the  main  techni¬ 
cal  sections.  A  more  complete  review  is  available  in  "Optics  of  the 
Atmosphere"  by  L.  ,1 .  McCartney  published  in  1976  (Reference  7).  W.  E.K. 
Middleton's  treatise  "Vision  Through  the  Atmosphere",  although  published 
in  1952,  is  still  a  good  introduction  to  the  subject  of  visibility  (Reference  8). 
Quantitative  discussions  of  visibility  will  be  given  in  Section  2.2. 

The  atmospheric  properties  that  are  important  to  visibility 
problems  are  those  governing  scattering  and  absorption  of  light  being 
propagated  through  the  atmosphere.  In  this  context,  the  atmosphere  is  an 
aerosol  which  is  defined  as  a  dispersed  system  of  small  particles  suspended 
in  a  gas.  The  gases  are  the  permanent  gases  such  as  nitrogen  and  oxygen 
and  the  variable  gases  such  as  ozone  and  water  vapor.  The  particles  come 
from  a  variety  of  natural  and  man-made  sources.  The  particle  size  range 
from  0.  01  to  10  microns  is  important  to  visibility.  Light  scattering  in  gases 
is  through  the  Rayleigh  mechanism  wherein  the  molecule  as  a  classical  dipole 
scatters  the  electromagnetic  wave.  Absorption  in  gases  is  by  the  molecular 
spectral  bands  which  are  generally  in  the  infrared  spectral  region.  Water 
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vapor  and  oxygen  have  near-infrared  bands  which  encroach  toward  the 

o 

visible  region  near  7000  A.  Water  vapor  bands  are  important  insofar 
as  they  are  to  be  avoided.  Oxygen  has  an  absorption  band  near  7  600  A 
called  the  A -band.  The  band  head  is  at  7  593  A  and  the  R  and  P  branches 
are  centered  around  7  619  A.  Another  weaker  band  called  the  B-band  is 

3  j 

centered  at  6900  A  with  band  head  at  6867  A.  In  a  later  section,  an 
attempt  is  made  to  use  the  A-band  for  discriminating  between  an  aerosol 
distribution  smoothly  decaying  with  height  and  a  distribution  which  has 
a  finite  layer  superposed  on  the  first  distribution.  Ozone  has  a  broad 
group  of  weak  bands  (Chappuis  bands)  between  5000  A  and  7000  A  which 
occasionally  need  to  be  considered  in  calculating  transmittance  for  optieal 
paths  passing  through  the  atmosphere.  They  are  not  further-  considered 
in  this  report. 

Rayleigh  scattering  of  electromagnetic  waves  by  molecules 
is  also  valid  for  small  particles  as  long  as  the  radius  is  much  smaller 
than  the  incident  wavelength.  (Near  the  peak  of  the  eve  sensitivity  curve, 
light,  has  a  wavelength  of  0.  55  microns.  )  As  the  particle  size  increases, 
the  interaction  between  wave  and  particle  becomes  very  complex  and  needs 
the  elaborate  theory  associated  with  G.  Mie  to  describe  the  scattering. 

Since  Mie  scattering  depends  on  particle  size,  scattering  quantities  for 
real  aerosols  are  obtained  by  averaging  over  particle  size  distributions. 
Absorption  in  Mie  scattering  is  described  phenomenologically  with  a  com¬ 
plex  index  of  refraction  whose-  imaginary  part  accounts  for  absorption. 

The  interaction  between  radiant  energy  in  the  form  of  electro¬ 
magnetic  waves  and  a  single  molecule  or  particle  is  quantitatively  described 
in  terms  of  a  cross  section.  The  product  of  this  area  and  the  power  per- 
unit  area  in  the  wavefront  is  the  amount  of  power  extracted  from  the  wave 

by  absorption  or  scattering.  To  account  for  a  large  number  of  molecules 

2 

or  particles,  this  cross  section  (dimension  L  )  is  multiplied  by  tin*  mole¬ 
cular  (or  particle)  density  (dimension  I,  S  to  produce  an  interaction 
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I, 


coefficient  with  the  dimension  of  an  inverse  length  ( L  ).  When  the 
interaction  is  absorption,  the  result  is  an  absorption  coefficient;  when 
the  interaction  is  scattering,  a  scattering  coefficient.  The  combination  of 
absorption  and  scattering  results  in  an  extinction  coefficient  which  is  the 
sum  of  the  two.  Those  scattering  coefficients  are  sometimes  further  identi¬ 
fied  as  being  a  Rayleigh  or  Mie  coefficient. 

These  symbols  have  to  be  annotated  by  subscripts  and  super¬ 
scripts  to  indicate  whether  the  interaction  is  absorption,  scattering,  or 
both  (i.e.  extinction)  and  to  identify  the  entity  (e.g.  molecules,  aerosols) 
which  perturbes  the  optical  beam.  To  avoid  notational  ambiguity,  it  is 
sometimes  necessary  to  display  what  appears  to  be  typographical  excess 
as  will  be  seen  in  Section  2.  2. 

Absorption  interactions  completely  remove  energy  from  the 
problem  and  can  be  described  solely  by  an  absorption  coefficient.  On  the 
other  hand,  after  a  scattering  event  energy  has  not  disappeared  but  has 
only  been  redistributed  in  direction.  (Redistribution  in  wavelength  can 
also  occur  but  is  not  important  to  this  exposition.  )  Since  this  redistribution 
is  seldom  isotropic  or  uniform  in  all  directions,  it  is  necessary  to  describe 
the  amount  scattered  into  the  various  directions  around  the  scattcrer.  This 
description  is  made  in  terms  of  the  amount  of  light  scattered  in  a  unit 
solid  angle  about  an  observing  direction  lying  at  an  angle  with  respect  to 
the  original  light  beam  direction.  This  quantity  can  be  expressed  cither 
as  an  angular  scattering  coefficient  or  as  an  angular  scattering  coefficient 
or  as  the  product  of  the  total  scattering  coefficient,  a  fraction  called  the 
phase  function  (the  term  phase  is  an  unfortunate  usage  which  is  analogous 
to  the  astronomical  use  of  phase  to  mean  angular  aspect,  as  in  the  phases 
of  the  moon),  and  a  normalization  constant  of  1/4  tt.  Using  the  same  symbols 
with  appropriate  subscripts  as  given  above  to  identify  the  specific  interaction 


mechanism  ,  uno  can  write  ' 


d  (60  =  .3  =  3  «M0)  (2.  1) 

4  n 

In  the  second  equality,  the  obvious  substitution  of  <t>  for  P/4  tt  is  made  to 
simplify  the1  notation.  The  name1  "phase  function"  will  also  be  used  for 
this  quantity. 

The  scattering  process  is  also  sensitive  to  the  polarization  of 
the  incident  electromagnetic  wave,  [Fnpolarized  light  (of  which  natural 
sunlight  is  an  example)  can  be  considered  to  be  the  sum  of  two  independent 
linearly  polarized  components  of  equal  intensity.  Thus  when  considering 
the  scattering  of  sunlight  on  molecules  and  particles,  the  scattering  co¬ 
efficient  or  cross  section  is  found  by  averaging  the  values  for  the  parallel 
and  perpendicular  components  of  the  beam. 

Rayleigh  scattering  is  described  by  a  well  known  formula  which 
can  be  used  to  calculate  scattering  coefficients  for  atmospheres  of  any  known 
composition.  Tabulations  of  phase  function,  total  scattering  coefficient, 
and  angular  coefficient  are  available  to  simplify  Rayleigh  scattering  com¬ 
putations  (see  particularly  References  7  and  9).  The  total  Rayleigh 

-27  2  -1 

scattering  cross  section  at  0.  55  microns  is  a  =  4.  6  x  10  cm  molec  .  Since 

ni  19  -3 

the  number  density  of  atmospheric  molecules  at  sea  level  is  2.  54  x  10  cm 

-7  -1  -3  -1 

the  total  scattering  coefficient  is  3  =1.17x10  cm  =  11.7  x  10  km 

m 

Values  at  other  wavelengths  are  readily  computed  using  the  well  known 
inverse  fourth  power  wavelength  dependence. 

The  calculation  of  Mie  scattering  by  a  single  particle  involves 
the  summation  of  a  slowly  converging  infinite  series  each  term  of  which 
is  a  complicated  function  which  must  be  evaluated  on  a  large  computer. 

For  a  mixture  of  particles,  these  calculations  must  be  repeated  many  times 
in  order  to  average  over  the  particle  size  distribution.  Needless  to  say, 

*  The  symbols  used  in  the  report  follow  that  of  Ref.  7  for  the  most  part.  The 
reader  is  cautioned  that  this  usage  is  not  necessarily  consistent  with  other  writers. 
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Particle  size  distributions  for  two  aerosol  models  (from  Reference  11). 
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Figure  2.  2.  The  phase  function  for  different  lower  atmospheric 
aerosol  models  at  a  wavelength  of  1.06  microns. 


NORMALIZED  EXTINCT 


re 


"A  ! 

10  L_ 
0.1 


•LOWTRAN  3 
CONTINENTAL 
MODEL 

- RURAL  AEROSOL  MODEL 


H 

-  H 


A 


-  URBAN  AEROSOL  MODEL 

135%  SOOT) 

- MARITIME  AEROSOL  MODEL 

{75%  SEA  SPRAY) 


TROPOSPHERIC  AEROSOL 


-1  L-LLLLlLL 


MODEL 

J _ i _ JJJJJil 


■I 


I  1 ..  J  1  1 1 1 1 


I.  10  100 

WAVELENGTH  (microns) 


2.  3.  Extinction  coefficients  for  several  aerosol 
models  (from  Reference  11). 


figure  2.4.  Vertical  distribution  of  attenuation  coefficient  at  0.  55/jm 

wavelength  for  several  aerosol  models  (from  Reference  11) 


The  presence  of  temperature  inversions  near  the  surface  will 
produce  a  region  of  nearly  constant  attenuation  coefficient  below  the  in¬ 
version  with  a  sharp  drop  above.  This  behavior  is  usually  associated  with 
low  values  of  meteorological  visibility  (e.  g.  2  to  10  km  as  shown  in  Figure 
2.4).  In  the  absence  of  inversions,  the  attenuation  coefficient  shows  a 
general  tendency  to  decrease  exponentially  with  a  scale  height  of  about 
1  km  as  shown  in  Figure  2.  5. 

The  way  that  these  properties  are  used  in  radiation  transfer 
problems  will  be  shown  in  the  next  section. 

2.  2  Radiative  Transfer 

Radiative  transfer  theory  describes  the  flow  of  radiation  through 
matter.  It  finds  extensive  use  in  that  part  of  astrophysics  dealing  with  the 
escape  of  radiation  from  stars.  The  calculation  of  the  brightness  of  tin*  sun 
illuminated  sky  is  made  using  the  equations  of  radiative  transfer.  The 
theory  that  applies  to  these  and  other  similar  problems  has  been  developed 
through  both  analytical  and  numerical  methods.  One  of  the  classic  expositions 
of  radiative  transfer  is  the  treatise  by  Chandrasekhar  (Reference  I'D  but  a 
more  readable  one  (but  not  less  mathematical)  is  that  by  Sobolev  (Reference 
14),  who  provides  much  material  on  the  transport  of  radiation  in  planetary 
atmospheres,  particularly  the  earth.  Not  covered  in  either  of  these  ref¬ 
erences  is  the  .Monte  Carlo  method  which  is  a  numerical  simulation  of  the 
entire  problem.  A  typical  Monte  Carlo  calculation  is  given  in  Reference  15, 
These  calculations  are  lengthy  and  therefore  do  not  replace  oilier  numerical 
calculations.  All  of  these  methods  of  solution  are  too  complicated  to  be  used 
to  make  simple  estimates  for  test  ing  concept s .  Some  recent  surveys  of 
methods  of  solving  radiative  transfer  problems  can  be  found  in  References 
16,  17,  and  18.  Fortunately,  Die  problems  treated  in  this  report  can  In- 
reduced  to  simple  forms  which  can  be  calculated  with  not  more  than  a  pro¬ 
grammable  hand  calculator. 


Relationships  of  Four  Aerosol  Scale  Heights  (H  ) 
With  Meteorological  Range  (V^  ),  Aerosol  Attenuation  P 
Coefficients  (p  )  and  Aerosol  Mixing  Layer  Altitudes  (0  -  5km). 
The  aerosol  scRle  height  family  was  computed  using  A  =  0.  55 id. 


The  dash  line  represents  values  of  0  (h)  above  5  km  (Elterman, 
1968)  P 


Figure  2.  5  Scheme  for  relating  ground  layer  aerosol  profile  with 
meteorological  range  (from  Reference  12). 


19 


The  physical  property  being  transported  or  propagated  through 
the  atmosphere  is  radiance,  frequently  called  by  the  imprecise  word  bright¬ 
ness.  This  quantity  is  the  amount  of  energy  emitted  in  unit  time  from  a  unit 

2 

area  of  a  source  into  a  unit  solid  angle.  The  cgs  units  are  watts  per  cm  per 
-2  -1 

steradian  (w  cm  sr  ).  (This  quantity  is  called  intensity  in  the  astrophysical 
literature  in  conflict  with  terminology  used  in  the  physics  literature.  )  The 
standard  symbol  used  for  radiance  is  N.  If  radiance  is  spectrally  depen¬ 
dent,  the  unit  wavelength  or  frequency  interval  is  added  to  the  units  and  the 
Greek  minuscules  A  (for  wavelength)  or  v  (for  frequency)  are  added  to  the 
symbol  either  as  a  subscript  or  in  parenthesis.  For  example,  if  the  spectral 

dependence  is  expressed  per  unit  Angstrom,  the  spectral  radiance  with  its 

-2  -1  -1 

units  would  be  written  N'f  A)  or  N  ,  (w  cm  sr  A  ).  Hadiance  is  the 

A 

appropriate  quantity  to  be  used  in  most  of  the  methods  to  be  discussed  be¬ 
cause  it  pertains  to  extended  sources  and  because  imaging  optical  systems 
preserve  radiance  (except  for  transmission  losses). 

Heuristically ,  the  equation  of  steady-state  radiation  transfer  is 
derived  by  considering  the  change  in  radiance  along  a  straight  path.  Changes 
occur  when  energy  is  lost  through  absorption  or  scattering  or  when  energy 
is  gained  through  emission  from  the-  atmosphere  or  from  scattering  of  energy 
into  the  direction  of  the  path.  Self-emission  of  the  atmosphere  is  onlv 
important  from  3  microns  and  beyond.  Thus  only  scattering  is  considered 
here. 

The  simplified  equation  of  transfer  can  be  expressed  in  words  as 

follows: 

f  radiance  change  along  a  unit  path]  =  |  radiance  decrease  along  unit  path 

from  absorption]  ■>-  |  radiance  decrease 
along  unit  path  from  scattering]  + 
[radiance  increase  along  unit  path  from 
scatte  ring] 
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bsing  tin-  atmospheric1  property  concepts  outlined  in  Section 
2.  1,  the  equation  of  transfer  ean  he  written  symbolically  in  the  form 


d  N 
d  x 
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T  (ft))  M (ft)  d  n  (2.  2) 


1'he  subscripted  coefficients  din  the  first  two  terms  on  the  right  hand  side 
of  the  equation  are  the  total  absorption  and  scattering  coefficients  for  the 
material  in  the  volume  element  around  the  path  ’length  dx.  Kuril  of  tile 
coefficients  can  be  broken  down  into  components  due  to  the  aerosol  and  the 
molecular  constituents.  1'hus 
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Because  of  the  linearity  of  the  absorption  and  scattering  loss  terms,  thesi 
coefficients  can  be  rearranged  to  form  extinction  coefficients  as  follows 
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Also  recognizing  that  scattering  light  from  molecules  is  by  the  Rayleigh 
mechanism,  the  subscripts  related  to  that  kind  of  scattering  ean  be  abbre¬ 
viated  to  "Rav".  Then  the  equation  of  transfer  can  be  rewritten 
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d  X 


In  many  useful  and  relatively  unrunijil  ieated  eases,  this  integro- 
differential  equation  ran  be  simplified  to  a  differential  equation  which  can  be 
easily  solved.  For  example,  bv  neglecting  the  in-scallering  term,  the  ex¬ 
ponential  transmission  formula  for  a  collimated  beam  can  be  immediately 
derived,  i.e. 


N  (x)  =  N  c 
o 


-(x-x  ) 
o 


(2.  4) 


where  N  is  the  beam  radiance  at  position  x  . 

o  o 

A  more  useful  example  can  be  obtained  by  considering  the 
apparent  radiance  of  an  object  having  an  intrinsic  radiance  N  located  at 
some  distance'  x  from  the  observing  point  in  a  sunlight  illuminated  atmos¬ 
phere.  On  making  the  reasonable  simplifications  in  liquation  2.2  that 
multiple  scattering  of  light  ran  be  neglected  and  that  an  effective  phase 
function  can  be  used  to  take  into  account  the  presence  of  Rayleigh  scat¬ 
tering,  the  following  simple  differential  equation  for  the  radiance  along  the 
ray  path  to  the  object  is  obtained. 


d  N  _ 
d  x 
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ext  seat  eft 


(2.  a) 


where  11  is  the  uniform  solar  irradianre  (W  cm  )  on  a  unit  path  length  on 
the  ray  path  and  d  is  the  scattering  angle.  The  solution  to  the  equation  is 
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(2.  G) 


When  x  o.  the  radiance  is  correctly  \  ,  and  when  x  -  ,  the  radiance  is 

o 

equal  to  (  .1  1  A  )  <t>  „„  II  which  is  Urn  radiance  of  the  horizon  seen  with 

ext  aer  eff 


with  no  intervening  object.  Thus  the  equation  can  also  lie  written 


N  =  N  (  1  -  e 
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(2.7) 


This  last  equation  is  used  to  determine  the  visibility  of  objects 
in  the  atmosphere.  An  object  is  visible  at  a  given  range  in  the  atmosphere 
if  the  difference  in  radiance  between  the  object  and  its  background  exceeds 
some  minimum  difference  appropriate  to  the  detector  being  used.  The 
value  of  the  distance  at  which  this  condition  holds  is  called  the  range  appro¬ 
priate  to  the  detector  under  consideration.  For  example,  when  the  eye  is 
the  detector  one  calls  the  distance  the  visual  range. 

If  the  contrast  of  an  object  is  defined  with  respect  to  the  horizon 
brightness  behind  the  object,  the  last  equation  can  be  manipulated  to  give 


N  -  N  “  (N 
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(2.  8) 


Dividing  both  sides  by  the  horizon  brightness  one  obtains  an 
equation  between  the  inherent  contrast  of  the  object  against  the  background 
and  its  apparent  contrast  when  viewed  at  a  distance  R 


C 


R 


C  e 
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(2.  9) 


It  is  useful  and  conventional  to  describe  the  clearness  of  the 
atmosphere  in  terms  of  a  distance  at  which  objects  are  no  longer  visible. 
The  most  common  such  definition  is  that  of  meteorological  range  which  is 
defined  as  that  distance  at  which  the  ratio  of  apparent  to  inherent  contrast 
at  a  wavelength  of  0.  55  microns  is  two  percent.  Thus  from  Kquation  2.  9 


r 


the  meteorological  range  becomes 

-In  (0.02)  3.912 

VR  =  3  =  ~  (2.  10) 

This  concept  of  contrast  and  contract  reduction  can  be  ex¬ 
tended  to  cover  more  general  situations  where  objects  are  viewed  with 
respect  to  backgrounds  other  than  the  horizon  and  where  contrast  is  relative 
to  another  object  instead  of  the  background.  An  application  of  the  extended 
concept  is  discussed  in  Section  4.  1. 
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In  the  context  of  this  report,  remote  sen  sing  is  tne  measure¬ 
ment  of  atrgpspheric  parameters  from  a  distance.  The  specific  concern 
is  with  measurement  of  atmospheric  visibility  and  its  component  parameters 
such  as  particle  size  distribution  and  particle  concentration  *  •  on.  ear**' 
satellites.  However,  it  is  best  to  recognize  at  the  outsei  tnat  «■  this  pc  ir.t 
in  time  it  only  appears  feasible  to  measure  the  total  vertical  transmission 
of  the  atmosphere  as  a  function  of  visible  wavelengths.  From  the  change 
of  transmission  with  wavelength  one  can  possibly  infer  ola--m-  01  particle 
size  distribution-.  Nor  a*-  this  time  does  it  seem  to  he  feasible  ro  mens  .re 
or  info-'  anything  about  the  variatior  of  any  parameter  with  jc  above 

■ :  ■  h  (■• . 

o .  ■  1 1  ’  i  i ;  i.  r  -  •  .■■■  i)e«m  am'  ming  us--d  to  is  ■  m.ner 
'-Pueiu  rnnpmt :  I'm  ■ -:o- •  -  ip :  r  Tea  m  r'  these  -,e 

i  .  *  ■  : '  pc  mam.  •.  e mi  c  !  te  ■  ; -r  r<  1<  w,v  •• 

■  por  no!  iJ  •  r  •  r  .  •-o'  i  ie  ■.  Ithoaru  ■  eLhou  s  ••  .  <  .  •<•  •> 

*■  - .  "•  tr  ■  .•  general)'-  m . - ■,  ltniue  t< •  the  ,u  .  j*v 

"*-’i  a  ’vi. . :*  or  hm  e  t,-.r  im  es  v.;T'  nrovi  i<  ->■;  ]  - 

or  •  i i c 1  i  •  •  Mir  •  i-’i'dlr.  v  r.i  'isi-ir  i : < -mot'-  -am.-  u_  torn:;  s< 

•  c.  ■■  r  m*'  dec  ‘ ,  rev  whir:  ;.o;- 

1  ''s.'..  v  m ■  ’  1 :  ■ r i  ■  e  rl  i  :  r  ■  h  a  s  l  r  T  e  r  a  U  t  e  1  v.  O  ■ :  t  a  ■  ■  .  t 1  - 

: r  •  ! :  •  ■  • :  ,  •  ••  -  ■  •  ’  •  .•  f  (  - •  •  •  .  •  .  .  i  - 


profile  is  calculated  by  inverting  an  integral  equation  which  connects  the 
measured  wavelength  dependent  radiance  values,  the  Planck  function  of 
wavelength  and  temperature,  and  a  weighting  function  varying  mainly 
with  altitude.  This  technique  has  been  applied  to  the  15  micron  IR  band 
of  CO  .  It  can  be  used  with  the  5  mm  microwave  band  of  O  This  same 
concept  has  been  used  to  obtain  vertical  profiles  of  water  vapor.  In 
principle,  the  technique  can  be  used  to  measure  the  ozone  vertical  pro¬ 
file  using  the  9.  6  micron  band;  in  practice  the  structure  of  the  band  is 
such  that  not  much  profile  information  can  be  obtained  although  the 
total  ozone  content  can  be  measured. 

Ozone  distribution  can  be  inferred  by  measuring  the  back- 
scattered  solar  (or  lunar)  UV  radiation  from  the  atmosphere  at  different 
wavelengths  in  UV  absorption  bands.  The  variation  in  absorption  changes 
the  effective  depth  to  which  the  radiation  penetrates  and  backscatters.  The 
changing  location  of  effective  scattering  allows  the  variation  in  ozone  con¬ 
centration  to  be  deduced. 

Measurements  on  aerosols  are  most  conveniently  carried  out 
in  the  wavelength  region  from  about  0.4  to  0.  9  micron.  In  this  range, 
there  is  little  absorption  or  emission  by  the  particles.  Thus,  the  possible 
techniques  which  can  be  used  to  measure  visibility  remotely  are  limited 
to  those  using  the  scattering  of  light.  A  simple  classification  scheme  for 
these  methods  is  illustrated  schematically  in  Figure  3.  1. 

The  elements  of  this  primitive  scheme  are  a  source  of  light, 
a  detector,  the  scattering  layer  whose  transmission  is  to  be  measured, 
and  a  boundary  which  may  scatter  light  diffusely  or  reflect  it  specularly. 

In  the  first  scheme,  the  source  and  detector  are  on  opposite  sides  of  the 
scattering  layer  and  the  transmission  is  measured  directly.  The  source 
and  detector  are  on  the  same  side  of  the  layer  in  the  second  scheme,  but 
the  detector  views  the  source  through  the  intermediary  of  a  specular 
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Figure  3.1.  Schematic  illustration  of  types  of  visibility  measuring 
techniques. 
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reflector.  In  the  third  scheme,  the  source  and  detector  are  again  on 
the  same  side,  but  the  intermediary  now  becomes  two  diffusely  re¬ 
flecting  surfaces  with  different  reflectances.  The  quantity  measured  is 
the  apparent  contrast.  The  detector  in  the  fourth  scheme  measures 
the  backscattered  light  from  the  layer;  the  diffusely  reflected  light 
from  the  boundary  surface  is  now  an  unwanted  background  term. 

No  actual  scheme  is  as  simple  as  these  skeleton  arrangements 
imply.  For  example,  most  actual  schemes  will  have  various  background 
or  otherwise  interfering  terms  for  which  corrections  must  be  made. 

Then  too,  the  scheme  must  be  elaborated  so  that  enough  independent 
measurements  are  made  to  allow  for  solving  for  the  unknown  variables 
of  the  problem.  For  example,  if  in  a  Type  1  experiment  the  source 
strength  S  is  known,  then  from  a  single  measurement  D  the  transmission 
of  the  layer  can  be  determined.  If  the  source  strength  is  unknown,  then 
a  second  measurement  is  required.  For  Type  1  experiments,  this 
additional  measurement  usually  takes  the  form  of  an  oblique  view 
of  the  source  through  an  effectively  greater  path.  If  the  scattering  layer 
is  a  thin  gas  the  second  transmission  is  equal  to  the  perpendicular  trans¬ 
mission  raised  to  the  power  sec  9,  where  9  is  the  angle  of  inclination  of 
the  second  path.  Algebraically  the  two  measurements  are  expressed 

Dj  =  T  Sj  and  D9=(T)SeC0S2  (3.1) 

Then  if  the  source  has  the  same  value  in  both  measurements,  the  trans¬ 
mission  of  the  layer  is 

r~~(r^2 

Vt)i 


^  1/ (sec0-l ) 


Rather  than  develop  similar  equations  for  the  general  types, 
the  various  specific  techniques  for  remotely  sensing  atmospheric  visibility 
will  now  be  considered. 
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PASSIVE  REMOTE  SENSING 


4.  1  Introduction 

This  section  presents  discussions  of  specific  methods  which 
have  been  proposed  in  the  literature  for  remotely  measuring  visibility 
from  satellites.  The  technical  principle  of  each  concept  is  presented. 

Also  included  in  the  discussions  are  the  advantages  and  disadvantages 
any  any  complicating  features.  Where  it  is  known  or  readily  estimated, 
the  sensitivity  is  discussed.  Any  results  or  pertinent  satellite  experiments 
are  extracted  from  the  original  presentation  and  presented. 

4.  2  Contrast  Reduction 

In  Section  2.2,  it  was  shown  that  contrast  reduction  along 
horizontal  viewing  paths  in  the  atmosphere  was  directly  related  to  the 
extinction  coefficient.  It  was  also  stated  that  the  concept  of  contrast 
reduction  could  be  extended  to  non -horizontal  paths  and  to  a  more  general 
background.  Such  an  extension  was  made  by  Duntley  (Reference  19)  and  is 
discussed  also  by  Middleton  (Reference  8).  Griggs  has  made  this  extended 
treatment  the  central  concept  in  a  method  for  measuring  the  total  aerosol 
content  of  the  atmosphere  (Reference  20). 

The  idea  of  using  a  change  in  contrast  of  objects  on  the  earth's 
surface  to  measure  haze  has  occurred  to  others,  for  example  O'Dell 
and  Weinman  (Reference  21),  They  predict  the  degradation  of  contrast  through 
the  use  of  numerical  calculations  of  optical  properties  of  a  sunlit  haze 
above  a  reflecting  surface.  In  their  approach,  O'Dell  and  Weinman  take 
the  reflectances  of  two  known  objects  along  with  the  mean  reflectance  over 
the  extended  surface  and  find  the  optical  thickness  of  the  haze  from  a  single 
equation  and  a  table  of  derived  factors.  They  also  point  out  that  if  the 
optical  depth  is  known  the  problem  can  be  inverted  and  unknown  surface 
reflectances  determined. 


Grigg's  approach  is  similar  in  that  ho  doveloped  a  model  to 
relate  optical  depth  of  the  atmosphere  and  the  apparent  contrast  of  the  two 
comparison  objects  whose  intrensic  contrast  is  known.  He  did  not  find  it 
necessary  to  account  separately  for  a  mean  reflectance  since  he  considered 
very  large  areas  (the  Salton  Sea  and  the  adjacent  desert  area).  His  model 
accounted  for  the  solar  zenith  angle  as  did  that  of  O'Dell  and  Weinman. 

Griggs  has  reported  determinations  of  aerosol  contact  from  ERTS-1  data 
taken  over  the  test  area  for  three  days  in  the  autumn  of  1972.  Lacking 
any  concomitant  ground  truth  data,  all  that  can  be  said  is  that  the  data  are 
reasonable . 

The  simplicity  of  the  method  and  the  (presumably)  ready 
availability  of  the  necessary  satellite  data  make  this  method  a  good  candi¬ 
date  for  further  exploration.  An  important  question  to  bo  answered  is 
how  accurately  must  the  inherent  contrast  of  the  comparison  objects  bo 
known. 

4.  3  Upwelling  Spectral  Radiance 

The  solar  radiance  reflected  from  the  earth  and  its  atmos¬ 
phere  depends  mainly  on  the  surface  reflectance  (or  albedo)  and  the  scat¬ 
tering  properties  of  the  atmosphere.  The  atmospherically  scattered 
radiation  depends  mainly  on  the  total  amount  of  scatterers  and  negligibly 
on  their  vertical  distribution.  Thus,  it  has  been  conjectured  that  in  cases 
of  low  surface  reflectance  it  would  bo  possible  to  deduce  the  amount  of 
aerosol  in  the  atmosphere  from  measurements  of  the  upwelling  spectral 
irradiance  from  the  earth -atmosphere  system  in  conjunction  with  a  numer¬ 
ical  model  of  thephysical  system.  This  approach  using  tin*  ocean  as  a  low 
reflectance  surface’  was  suggested  by  Griggs  (References  20  &  22)  who  initially- 
reported  tentative  tests  using  ERTS  satellite  data.  Subsequentially,  he  has 
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reported  further  observations  using  data  from  the  Landsat  series  satellites 
(Reference  23)  and  from  the  NOAA and  GOES  satellites  (Reference  24).  All  of 
his  results  show  over  oceans  there  is  a  linear  relation  between  the  upwelling 
visible  radiance  and  the  atmospheric  aerosol  content. 

A  more  extensive  test  of  the  concept  has  been  reported  by 
Koepke  and  Quenzel  who  used  data  from  the  geostationary  satellite  SMS 
taken  at  the  time  of  Project  GATE  (Reference  25).  The  model  radiance  was 
calculated  with  a  numerical  method  including  multiple  scattering  and  a 
realistic  model  of  the  ocean  surface  including  wave  structure.  Figure  4.  1 
shows  their  contour  plot  of  atmospheric  turbidity  in  the  Atlantic  Ocean  from 
the  Caribbean  to  the  African  Coast  on  31  July  1974.  The  plot  also  includes 
ground-truth  measurements  from  stations  operated  for  Project  GATE. 

Figure  4.  2  shows  the  variation  of  measured  radiance  at  a  fixed  location 
(Dakar)  throughout  the  day.  Superimposed  on  the  satellite  data  are  calcu¬ 
lated  radiance  [or  two  aerosol  models,  one  of  which  includes  Sahara  dust. 

While  this  method  has  been  demonstrated  to  work  over  oceans 
it  does  not  seem  to  have  been  explored  over  land  masses.  The  next  logical 
step  would  be  to  find  out  whether  it  works  over  soil  and  vegetation  regions 
of  low  reflectance.  It  has  been  suggested  that  some  continental  lakes  could 
be  used  as  low  reflectance  backgrounds  (Reference  26). 

4.4  Simulated  Polar  Nephelometer 

An  interesting  method  for  obtaining  the  optical  thickness  of  the 
earth's  atmosphere  from  staellite  based  observations  has  been  proposed 
by  Livshits  et  al  in  Reference  28.  The  principle  is  similar  to  that  used 
to  obtain  scattering  coefficients  from  polar  nephelometer  data.  If  the 
scattering  distribution  function  g  ( <j>)  =  (3<P(<f>)  is  known,  then  the  scattering 
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To  obtain  a  working  equation  which  can  be  used  to  process  satellite  obser¬ 
vations,  the  following  argument  is  made. 

The  observed  scattered  light  intensity  has  components  due  to 
first-order  and  multiple  scattering  and  to  ground  reflected  light.  Consider 
the  difference  of  two  intensities  measured  symmetrically  about  the  zenith 
in  a  vertical  plane  through  the  sun.  If  the  reflectance  of  the  surface  is 
symmetrical  about  the  zenith  (a  weaker  condition  than  requiring  a  Lambertian 
surface),  the  intensity  difference  contains  mainly  the  first-order  scattering 
component.  An  additional  assumption  which  must  be  made  is  that  back- 
scattered  light  at  scattering  angles  greater  than  90°  can  be  neglected.  If 
these  assumptions  hold,  then  the  intensity  difference  for  a  given  zenith  angle 
is  approximately  equal  to  the  aerosol  scattering  function  R  I*  (^)  . 
Measurements  covering  all  scattering  angles  can  then  be  used  to  find  the 
aerosol  optical  thickness  by  integrating  over  the  scattering  angles.  In 
practice,  it  is  difficult  to  find  the  right  conditions  to  make  the  necessary 
satellite  observations.  In  particular,  the  important  forward  scatter  contri¬ 
bution  will  be  hard  to  get. 

The  working  equation  given  in  Reference  28  (symbols  have  been 
modified  to  agree  with  our  notation)  is 

(  t  +  -  )  (m.  +  m  ) 

ner _ R _ 1 _ s _ 

<-(  r  +  r  )  (m.  +  m  )  )|  m. 
aer  R  1  s  J  i 

sin  <t> .  (4.1) 
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The  diagram  in  Figure  4.  3  is  useful  in  understanding  the  notation.  N  ^  and 
N'9  .  are  the  radiances  of  the  earth's  atmosphere  observed  from  two  directions 
making  angles  0  .  symmetrically  located  about  the  zenith.  The  plane  con¬ 
taining  these  two  directions  lies  in  the  plane  formed  by  the  zenith  direction 

and  the  sun  direction.  The  scattering  angle#  is  defined  by  the  direction  of 

i 

radiance  N  and  the  sun  direction.  The  aerosol  vertical  optical  thickness 


of  the  atmosphere  is  t  and  the  Rayleigh  optical  thickness  *  .  The  muss 
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of  the  atmosphere  in  the  direction  of  the  sun  is  mg  and  in  the  direction  of 

the  two  radiances  m;  .  II  is  the  solar  illuminance  outside  the  atmosphere 
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and  perpendicular  to  the  sun  direction,  a  «J»  is  the  difference  in  the  Rayleigh 
phase  function  at  the  two  observation  directions.  In  the  form  given  here,  the 
integral  is  approximated  by  a  discrete  sum  with  n  points  and  an  angular  inte¬ 
gration  step  h  in  radians.  An  iterative  p-ocess  must  be  used  to  find  Taer 
from  the  given  observation  of  radiance  in  units  of  the  solar  illuminance  H 
and  calculated  values  of  t  and  the  phase  function  <J>^. 

The  derivation  of  the  working  equation  was  not  given  in  Reference 
28  which  was  only  a  brief  note.  Since  the  source  of  the  equation  was  given 
as  an  untranslated  Russian  reference  not  readily  available,  an  independent 
derivation  is  given  in  Appendix  A.  This  derivation  disclosed  two  other  im¬ 
plicit  assumptions:  (1)  that  r  and  have  the  same  vertical  altitude 
dependence,  and  (2)  that  the  phase  functions  <l>  and  <J>  are  constant  with 
altitude.  Neither  of  these  are  true  in  real  atmospheres,  but  their  effect  may 
not  be  significant  in  practice. 

This  concept  can  also  be  applied  to  the  inverse  situation  where 
observations  are  made  from  the  ground.  A  working  equation  for  this  situation 
was  also  given  in  Reference  27  and  appears  in  Appendix  A  as  Equation  A.  18. 
Optical  thicknesses  obtained  from  this  ground  based  method  were  used  to  cal¬ 
culate  transparency  coefficients  (vertical  atmospheric  transmission)  which 
were  compared  to  transparency  coefficients  obtained  by  the  Bouger  method 
(Langley  method).  The  correlation  between  the  two  results  was  good  as  can 
be  seen  from  Figure  4.  4  which  is  reconstructed  from  the  original.  (The 
optical  wavelength  used  was  not  specified).  Although  it  is  not  explicitly  stated, 
it  is  apparent  that  these  comparisons  were  used  to  find  a  correction  factor  to 
account  for  missing  angular  data. 

3fi 


Data  from  the  Kosmos  149  satellite  were  used  to  test  the  method 
given  by  Equation  4.  1.  Only  those  data  which  satisfied  the  following  conditions 
were  used.  First,  the  solar  zenith  distance  had  to  be  large  enough  that  an 
adequate  range  of  scattering  angles  was  obtained  yet  not  so  large  that  the 
earth's  curvature  effects  were  significant.  A  solar  zenith  distance  of  75° 
was  found  to  be  suitable.  Second,  the  measurements  had  to  be  in  (or  near) 
the  plane  of  the  sun's  vertical.  Finally,  only  data  from  a  cloudless  atmosphere 
were  valid. 

Data  satisfying  all  of  these  criteria  were  used  with  the  working 
equation  to  find  a  value  of  optical  thickness  which  was  then  corrected  by  an 
empirical  correction  factor  that  had  been  previously  determined  by  comparisons 
between  other  values  determined  by  this  method  and  ground  measurements  of 
the  Bouger  optical  thickness  (vertical  atmospheric  thickness  determined  by 
Langley  method).  No  ground  truth  data  was  available  to  validate  these  results. 
The  authors  state  that  the  histogram  of  optical  thicknesses  obtained  from  the 
satellite  technique  was  similar  to  those  obtained  in  ground  measurements. 

They  further  add  that  the  results  suggest  that  the  asymmetry  of  the  reflected 
light  had  little  influence  on  the  determination,  in  agreement  with  the  assumption. 

In  spite  of  the  several  assumptions  whose  validity  can  be  questioned, 
this  method  is  potentially  important  in  that  it  removes  the  ground  reflectance 
from  the  problem.  At  the  same  time  the  technique  of  subtracting  two 
measurements  is  its  main  weakness,  since  the  residual  difference  must  be 
greater  than  noise  and  instrumental  fluctuations  in  the  two  symmetrical 
measurements.  In  practice,  this  restriction  will  limit  measurements  to 
regions  of  low  but  finite  reflectance. 
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4.  5  Nighttime  Observations  of  Cities* 


The  vertical  optical  thickness  of  the  earth's  atmosphere  is  fre¬ 
quently  determined  from  ground-based  stations  by  a  technique  called  the 
Langley  method  in  American  literature  and  the  Bouger  method  in  European 
literature,  A  photometer  is  aimed  at  the  sun  and  its  output  is  recorded  as 
a  function  of  the  sun's  zenith  angle.  Using  the  approximation  that  the  path 
through  the  atmosphere  is  proportional  to  the  secant  of  the  zenith  angle  for 
angles  less  than  80  degrees,  the  output  signal  is 
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(4.  2) 


and  the  natural  logarithm  of  the  signal  is 


In  S  =  In  S  -  T  sec  0  (4.  3) 

o  o 

or  the  common  logarithm 


log  S  =  log  Sq  -  sec  0  log  e  (4.  4) 

=  logS  -  0.4343  t  sec  0 
o  o 

A  plot  of  the  common  log  of  the  output  against  sec  0  is  a  straight  line  whose 
intercept  is  log  and  slope  is  (-0.  4343  t  ).  The  vertical  thickness  of  the 
atmosphere  r can  be  calculated  from  the  slope. 

This  method  could  be  reversed  if  there  were  isotropic  light 
sources  on  the  ground.  The  illumination  of  modern  cities  at  night  is  intense 
enough  to  be  detected  but  is  most  certain  not  to  be  isotropic.  Excellent  photo¬ 
graphs  of  the  nighttime  earth  can  be  found  in  References  29  and  30. 


*T he  concept  in  this  section  developed  during  discussions  with  E.Shettle 
of  an  earlier  idea. 


The  angular  variation  of  the  un-re  solved  light  source  can  be  resolved  into 
Fourier  components  which  will  not  be  known  a  priori  but  might  possibly  be 
determined  from  the  data. 


To  examine  the  feasibility  of  this  concept,  a  hypothetical  source 
composed  of  an  isotropic  source  and  a  cosine  varying  source  in  different 
proportions  will  be  investigated.  The  intensity  of  the  source  as  observed 
from  outside  the  earth's  atmosphere  as  a  function  of  zenith  angle  will  be 
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By  normalizing  to  the  value  at  the  vertical  and  writing  cos  0  -  u,  the 
expression  becomes 
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Further  mathematical  manipulation  can  put  this  relation  in  a  more  convenient 

form.  Let  d(\+c)  d  and  ll(\+c)  -  1  -  d.  Also  let  .v  (—  -  11  -  (1-ul/u 

U 

Then  ;(  r  w)  -  f  i  -  d( - r)  ]  e  .  The  constant  d  measures  the  nro- 

o  w  *  1 

portion  of  cosine  variation  in  the  source  and  ^anges  from  zero  ta  unitv.  When 

the  zenith  angle  varies  from  0°  to  85°,  the  variable  w  ranees  from  0  to  10. 

Figures  4.  5  to  4.  7  show  how  the  normalized  intensity  varies  as  a  function  of 

angle  through  the  variable  w  for  different  values  of  the  vertical  optical  thickness 

r  ,  The  parameter  d  varies  from  curve  to  curve,  <ince  examination  of  these 
o 

graphs  does  not  reveal  anv  obviously  simnle  wav  of  recuverine  the  two  unknown 
quantifies  -  and  n,  rim  erica!  methods  must  lie  used. 
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The  working  equation  giving  the  measured  intensity  as  a 

function  of  w,  r  and  d  can  be  alternately  solved  for  7  and  d 
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If  the  intensity  is  measured  at  two  different  angular  positions,  a  simple 

iterative  numerical  process  (described  in  Appendix  B)  can  be  used  to  solve 

for  r  and  d.  This  procedure  has  been  used  to  estimate  the  precision  to  which 
o 

these  values  can  be  determined.  The  exact  values  of  j  are  systematically 

varied  by  t  1  percent,  and  these  perturbed  values  then  used  to  calculate  values 

of  •  and  d.  Table  4.  1  shows  the  resulting  variations  in  these  values.  Since 
o 

the  variations  appear  to  be  linear,  assumed  variations  of  _  T>  percent  in  the 

measured  intensities  would  cause  variations  in  the  calculated  r  and  d  to  be 

o 

five  times  as  great.  Thus,  in  the  worst  case,  there  will  be  an  error  of  about 
30  percent  in  the  value  of  vertical  optical  thickness  deduced  from  a  pair  of 
measurements.  This  error  can  be  reduced  by  averaging  the  results  of  many 
pairs  of  measurements.  Although  this  analysis  needs  to  be  extended  to  addi¬ 
tional  terms  in  the  Fourier  expansion,  it  shows  considerable  promise  at  this 
point. 

It  is  probable  that  the  intensity  and  angular  distribution  of  light 
from  a  given  city  remains  approximately  the  same  from  night  to  night.  In 
that  case,  the  city  would  have  characteristic  values  of  "j  "  and  "d"  for  each 
hour  of  the  night  (excluding  perhaps  deviations  resulting  from  night  baseball 
or  football  games,  major  fires,  etc.  ).  Such  characteristic  values  could  be 
used  as  initial  guesses  in  the  iterative  solution  or  even  as  the  basis  for 
estimating  the  optical  depth  from  a  single  measurement. 
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4.6.  Reciprocal  Observations  of  Differences  in  Brightness 

In  section  4.2,  a  method  for  obtaining  the  transmission 
of  the  atmosphere  by  measuring  the  reduction  in  contrast  of  two  areas 
on  the  earth  with  different  reflectivity  was  described.  It  is  a  Type  3 
method  in  our  simple  classification  scheme  and  requires  a  single 
observation  along  with  a  prior  knowledge  of  the  two  reflectances.  This 
method  is  generalized  in  this  section  to  get  a  method  which  does  not  need 
values  of  reflectances,  but  at  the  price  of  complexity  in  number  and 
condition  of  observations. 


Il 


A  surface,  a,  is  irradiated  by  attenuated  direct  sunlight 
and  forward  scattered  sunlight.  The  direct  term  is 


m  ~  T  sec  0 

H  cos  <f>  ,  e  1 

o  1 


(4.9) 


where  0j  is  the  zenith  angle  of  the  sun  and  t  is  the  optical  depth  of  the 
atmosphere  evaluated  along  the  vertical. 

The  forward  scattered  radiation  reaching  the  surface  is 

-  t  sec  <t> 
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where  (  1-e  ) 


is  the  total  of  scattered  light  and  R  (  0  )  is  the 


fraction  of  that  total  which  is  scattered  in  a  downward  direction. 

The  inherent  radiance  of  the  surface  N  (0,0)  when  viewed 

a  1  ? 

at  angle  0  (measured  from  the  zenith)  is 
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where  is  the  reflectivity  of  surfaces  and  rj^f  (0^0^)  is  the  probability 

that  a  photon  incident  at  angle  0  will  be  scattered  into  an  element  of 

1  -1 

solid  angle  in  the  direction  0  (sr  )  and  p  g  (ft  0  )is  the  probability 

Z  3  &  1  Z 

that  forward  scattered  radiation  incident  on  the  surface  will  be  scattered 
into  an  element  of  solid  angle  in  the  direction  09  (sr'V 
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The  apparent  radiance  of  the  surface  N  i4>  <f>t)  when  viewed 

u  12 

from  outside  the  atmosphere  is 


N  (<*><)>  )  =  N'  e 

«  1  Z 
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where  U  (  is  the  upwelling  light  which  has  been  scattered  in  the  atmos¬ 

phere. 


The  apparent  radiance  of  a  second  surface,  b,  adjacent  to  a  may 
be  described  as  follows: 

- t  sec  0 
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U-e  *)  R  (0:)  1  +  U  (01  02) 


The  brightness  difference  between  surfaces  a  and  b  when  illuminated  from 
and  viewed  from  <f> 2  is  where. 


^b^lty  =  Hocos01|e 
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The  brightness  difference  between  surfaces  a  and  b  when  illumi¬ 
nated  from  <t>2  and  viewed  from  <f>  ^  is  D_u(  where. 


Dab(<^2<?>l)  =  Hocos02fe 
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The  coefficients  and  f^  (</>^</>2^  as  use<^  here  meet  the  requirements 


(4.  16) 
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of  the  reciprocity  theorem  and  hence 
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So  far  no  assumptions  have  been  made.  To  go  further  requires  some 
relationships  between  R{  <£j)  and  R{^2)<  between  g&(  9)  a.nd 

and  between  g^(  ^  and  g^(il  2  ^  1 )  . 

The  coefficient  B  (<f>  )  has  to  do  with  the  scattering  of 
sunlight  incident  from  zenith  angle  <t>  .  R  (<f>)  is  the  fraction  of  the 
scattered  sunlight  which  continues  in  a  downward  direction.  This 
fraction  is  a  function  of  the  zenith  angle  <£,  If  a  model  for  the  scattering 
process  is  assumed,  for  instance,  Rayleigh  scattering  plus  an  appropriate 
quantity  of  a  continental  aerosol  as  described  by  Deirmendjian  then,  for 
any  value  of  r  ,  the  optical  depth  of  the  vertical  atmosphere,  the  quantity 
of  aerosol  is  fixed  and  R  (4>)  can  be  calculated  easily  for  all  values  of  (4>) 
as  long  as  single  scattering  dominates.  In  cases  for  which  aerosol 
scattering  is  reasonably  greater  than  molecular  scattering,  which  is 
always  the  case  for  the  red  end  of  the  visible  spectrum  and  usually 
the  case  for  the  total  visible  spectrum  except  in  very  clear  weather, 
the  ratio  R  (4>)  can  be  fixed  for  any  value  of  <f>  without  regard  to  r  . 

At  this  stage  it  is  assumed  that  this  is  the  case  and  values  of  R  (<t>)  can 
be  assigned.  A  first  cut  at  calculating  R  (<t>)  is  presented  in  Appendix  C. 
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The  coefficients  such  as  $2^  havc  to  <4°  with  4he  probability 

that  a  photon  incident  on  surface  a  from  a  direc  tion  within  the  solid 
angle  fi  j  with  central  axis  in  the  direction  will  be  scattered  into  an 
eleme  nt  of  solid  angle  in  the  direction  <^>2  •  The  distribution  about  0^ 
of  the  numbers  of  photons  within  n  ^  will  be  quite  similar  to  the  phase 
function  for  scattering  and  will  have  a  strong  forward  lobe.  With  this 
in  mind  one  can  write 


ga<V2> 
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since  the  conditions  for  reciprocity  are  approached  but  not  met.  If 
the  surface  a  is  Lambertian  then. 
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where-  pa  is  the  reflectivity  of  the  surface  a  .  In  what  follows  it  will 

be  assumed  that  g  ( Q  <j>  )  g  (n  0  )  for  whatever  reason.  At  this 
a  12  ‘a  21 

stage  equation  4.21  will  be  used. 

Equation  4.  19  can  now  be  rewritten 
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Now  imagine  that  a  second  pair  of  observations  are  made  involving  angle 
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and  finally 
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One  notes  that  the  left  hand  side  of  equation  4.  24  involves  measurable 
quantities  and  the  right  hand  side  involves  only  one  unknown,  r  ,  In  prin¬ 
ciple,  then,  if  the  assumptions  made  are  acceptable,  the  optical  depth  of 
the  vertical  atmosphere  can  be  calculated  from  four  measurements  of  the 
difference  in  brightness  of  two  adjacent  areas  on  the  earth's  surface.  These 
measurements  are  not  made  casually.  Exacting  relationships  to  the  solar 
zenith  angles  are  defined  by  the  conditions  discussed  in  the  text  .hove. 


A  variation  on  the  method  is  of  some  interest.  Imagine  that 


the  surface  a  is  of  great  extent  and  the  surface  b  is  small  and  a  ell  within  the 
boundaries  of  a.  For  example,  a  might  extend  beyond  b  at  least  10  times  h, 
the  scale  height  of  the  aerosols  (  h  «  1  km)  and  b,  while  resolvable,  might 
have  no  dimension  greater  than  0.  2h.  Under  these  conditions  b  is  a  small 
perturbation  in  a.  If  in  the  atmosphere  above  the  surface  all  gradients  in 
scattering  effects  are  vertical  (an  assumption  made  throughout  this  section) 
then  for  regions  of  a  a  few  kilometers  in  from  its  edges  the  reciprocity  principl 
will  hold  for  the  combination  surface  plus  atmosphere.  That  is 
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The  left  hand  side  of  equation  4.  24  may  be  written 
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using  equations  4.  25  this  becomes 
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This,  now,  may  be  set  equal  to  the  right  hand  side  of  equation  4.24.  In 
this  interesting  case  only  the  brightness  of  the  small  patch  b  must  be 
measured  from  four  positions  to  obtain  r  .  Again  these  positions  and 
their  relationships  to  the  zenith  angle  of  the  sun  are  established  in  a 
rather  exacting  way  by  the  foregoing  text. 


The  observing  position  requirements  cannot  be  satisfied  by  a 
single  orbiting  satellite  progressing  in  same  direction  as  the  sun.  While 
it  may  not  be  cost  effective  to  devote  multiple  satellites  to  routine  obser¬ 
vation  of  visibility,  this  method  might  be  employed  for  one-time  calibration 
or  comparison  purposes.  It  could  also  be  used  in  those  fortuitous  cases 
where  the  observing  positions  occur  with  existing  satellites  in  orbit. 
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4.7  Upwelling  Spectral  Radiance  Observed  in  Oxygen  A  Band 

A  major  source  of  difficulty  in  attempting  to  measure  the  amount 
of  atmospheric  aerosol  from  the  upwelling  spectral  radiance  as  discussed 
in  Section  4.  3  is  the  presence  of  light  reflected  from  the  earth's  surface. 
The  methods  given  in  Section  4.  3  were  limited  to  conditions  where  this 
component  of  the  measured  light  was  small  and  was  accounted  for  through 
a  model  of  the  process.  The  method  of  Section  4.4  eliminates  the  ground 
reflection  by  subtracting  two  signals,  one  of  which  is  essentially  the  back¬ 
ground  reflection.  The  possibility  of  eliminating  the  background  by  absorp¬ 
tion  is  explored  in  this  Section.  Since  the  proposed  method  does  not 
appear  to  have  sufficient  sensitivity,  only  an  outline  of  the  concept  and 
the  final  results  are  presented. 

It  was  mentioned  in  Section  2.  1,  that  oxygen  molecules  have  an 
absorption  band  with  origin  at  0.  7619  microns.  If  the  earth  is  observed 
in  narrow  spectral  resolution  through  this  band,  then  different  amounts 
of  total  path  absorption  will  occur  as  different  parts  of  the  band  are  used. 
And  since  the  surface  and  the  aerosol  scatterers  are  located  at  different 
places  on  the  light  path,  the  contributions  of  each  to  the  reflected  light 
should  be  weighted  differently.  In  principle  then,  it  might  be  possible  to 
separate  out  the  surface  reflectance. 

To  investigate  this  concept,  the  following  simple  problem  was 
set  up  and  evaluated  on  a  HP-19C  programmable  calculator.  Radiation 
from  the  sun  is  incident  on  the  earth's  atmosphere  and  surface  at  a  given 
zenith  angle.  The  single  scattered  reflected  spectral  radiance  is  calculated 
at  various  narrow  intervals  in  the  oxygen  A-band.  It  is  further  assumed 
that  strong  forward  scattering  effectively  eliminates  any  scattering  ex¬ 
tinction  effects  for  both  the  incoming  solar  flux  and  the  upwelling  single 


scattered  flux.  The  Elterman  model  for  different  surface  meteorological 
ranges  (Reference  12  and  Figure  2.  5)  was  used  to  describe  the  atmospheric 
aerosol.  Calculations  were  made  for  one  surface  reflectance  (  15  percent) 
two  viewing  configurations  and  two  surface  meteorological  ranges  (5  km  and 
20  km). 

Figure  4.  8  shows  the  results  with  sun  at  45°  zenith  angle  and 

the  observer  at  the  zenith,  while  Figure  4.  9  shows  results  with  the  sun 

again  at  45°  but  with  the  observer  now  at  90°  with  respect  to  the  sun.  Each 

graph  shows  the  total  radiance  N  and  its  components  from  the  aerosol,  N  : 

aer 

Rayleigh  scattering,  N  ;  and  the  ground  N  .  The  range  of  absorption  co- 

R  G 

efficient  is  sufficient  to  reach  the  bottom  of  the  absorption  band.  Although 
there  is  a  calculational  difference  between  the  two  visibility  conditions,  it 
does  not  appear  to  be  large  enough  to  distinguish  between  them. 

It  is  possible  that  over  oceans  when  the  surface  reflectance  is 
small,  then  some  vertical  resolution  of  the  atmospheric  attenuation  may  be 
achieved  by  this  method.  However,  over  land  where  surface  albedos  range 
from  0.  1  to  0.  5  for  surface  areas  and  up  to  0.  9  over  clear  snow,  the 
expectation  of  achieving  any  vertical  resolution  is  not  high.  Further  com¬ 
plications  arise  because  the  0.7  6  micron  band  falls  within  the  high  reflecti¬ 
vity  chlorophyll  band  and  strong  seasonal  variations  in  surface  reflectivity 
can  be  expected. 

Subsequent  to  our  analysis,  Badayev  and  Malkevich  published  a 
paper  describing  a  similar  concept  (Reference  34)  but  using  a  more  sophis¬ 
ticated  mathematical  approach.  They  consider  the  method  to  be  useful  over 
weakly  reflecting  surfaces  on  the  basis  of  their  numerical  calculations. 

Also  they  appear  to  be  more  optimistic  about  making  corrections  for  multiple 
scattering  and  finite  surface  reflectances. 


Absorption  rneffuient ,  a  km' 


Figure  4.8.  Calculated  spectral  radiance  of  the  atmosphere  in  different 
parts  of  the  oxygen  A -band  (characterized  by  absorption  coefficient  a) 
seen  from  above  for  the  conditions  given  on  the  graph. 
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5. 


ACTIVE  REMOTE  SENSING 


The  previous  section  has  described  ways  of  measuring 
the  aerosol  content  of  the  lower  troposphere  using  mainly  scattered 
sunlight  but  am  '  considering  the  use  of  Light  sources  on  <-amh  at 
night.  Anotner  possible  way  of  measuring  these  aerosols  is  to  use 
active  light  sources  in  artificial  satellites.  The  only  feasible  source 
is  a  pulse  laser  used  in  a  iidar  configuration.  Both  technical  and  political 
considerations  affect  the  use  of  this  source.  The  important  technical 
factors  are  theoretical  feasibility,  engineering  practicality,  and  safety. 
Lidars  nave  been  in  use  for  many  years  for  measurement  of  aerosols 
along  paths  near  the  surface  of  the  earth  and  for  measuring  the  vertical 
distributions  of  aerosols  along  vertical  paths  extending  to  the  stratospnere 
The  theoretical  principles  of  their  use  have  been  weLi  established.  Two 
studies  of  the  feasibility  of  using  satellite  based  lidars  to  measure  the 
content  of  aerosols  and  other  constituents  in  the  atmosphere  are  given 
m  Reterences  31  and  32.  While  Reference  31  limits  its  consideration  to 
aerosoLs  above  10  km  and  does  not  address  our  prooiem,  it  does  provide 
much  informative  review  and  discussion. 

Reference  :  1  presents  a  detailed  discussion  of  lidar  system 
whi  m  might  be  useful  to  measure  aerosol  distribution  from  the  proposed 
Space  Shuttle.  They  conclude  that  an  eye  safe  svstern  operating  from  an 
altitude  of  185  km  is  theoretically  possible  but  that  it  was  at  the  Emil 
ot  i-u  rr  •  nt  teehm  >7oev  ( 15*7 "A.  The  svstern  operates  in  high  spectral 
resolution  (0.002  A)  both  at.  the  transmitter  and  roeer'  r .  Fabry -P'  rot 
eta  ions  ire  usee  to  obtain  this  resolution.  The  source  was  a  N,,  sun.  ;a- 
<iiam  VE1  A)  ia  Dumping  a  dvo  (4000 -5000  A'  in  a  laser -etalon-anpli- 
m'jgu  v-i '  i‘  • '.  Th"  ern  rgv  per  Dulse  was  50  u  i  .  '  •  p  iis. 

’  ■.  an  1  'h  rept  tit  ion  rate  was  no  TEt. 
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It  is  of  interest  to  compare  this  proposed  system  with 
the  system  that  is  currently  being  suggested  for  Space  Shuttle  use 
(Reference  33).  Table  5.  1  taken  from  this  reference  summarizes  para¬ 
meters  needed  to  trace  global  flow  of  pollutants  including  aerosols. 
There  are  to  be  other  measurements  of  other  parameters  reu*ic~.  ^ 
meteorology,  chemical  composition,  transport  of  thermospheric  atomic 
species  and  others.  A  standard  hardware  set  to  be  used  for  all  experi¬ 
ments  is  given  in  Table  5.  2.  Finally,  the  outline  of  the  experiment  is 
given  in  Table  5.  3.  It  will  be  of  great  interest  to  follow  the  progress  of 
these  proposed  active  remote -sensing  experiments. 


Table  5. 1.  Parameter  Summary  Extracted  from  Reference  33 


SUMMARY  OF  PARAMETERS  TO  BE  MEASURED  TO  TRACE 

global  flow  of  pollutants 


parameter 

Height 

Range 

Ax 

Ax 

Aerosols 

0-10 

100 

1 

10-30 

500 

2 

h2o 

0-10 

100 

1 

10-20 

500 

2 

total 

100 

°3 

0-10 

100 

1 

10-30 

500 

2 

XYre.g.  ,S02 

0-10 

100 

1 

gh4  nh3,no2. 

10-50 

500 

3 

N„0,  CFMs,  etc. 

it 

total 

100 

Temperature 

0-3 

100 

0.  2 

5-20 

500 

1 

\J 

0-10 

100 

1 

V 

X 

10-50 

500 

2 

0-10 

100 

1 

Vz 

10-50 

500 

2 

CO„  mixing  ratio 

0-3 

100 

3 

Remarks 

20%,  haze  layers, thin 
cirrus, Saharan  &  Mon¬ 
golian  dust 
20%,  stratospheric 
aerosols 

20  % 

20  % 

20  %_ 

20  % 

20  °fo 

20  %' 

20  % 

20  % 

IK,  Define  boundary 
layer 

2K,  Define  tropopause^ 

5  ml  s  j 
10  m  /  s  J 

5  cm/s 
10  cm/s 

1%,  Identify  sources 
and  sinks 


Exp.  No- 
2,3 

6,  24 
9,10 

10,  12 
18,  17 

10,  18,23 

17 

19.  20 
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Table  5.2.  Lidar  System  Parameters  Extracted  from  Reference  33 


ASSUMED  STANDARD  HARDWARE  SET  FOR  EXPERIMENT 

DESCRIPTIONS 


Lasers: 

Nd-YAG  1  .1/pulse;  10  pps,  TEM  00  mode;  15  nsec  pulse; 

^  1  A  linewidth;  5%  amplitude  stability 

Derivatives  X2  -  35  % ;  X4  -  10  % 

600  nm  (dye)  -  10%  (0.05  A)  or  30%  (1A) 

300  nm  (dye  X2)  -  1%  (0.  05  A)  or  3%  (1A) 


^^2  1  M  pulse;  300  pps;  TEM  00  mode;  100  nsec  pulse; 

10  MHz  stability  and  linewidth  in  9-12  pm  range  on 
selected  lines;  continuous  turning  and  other  lines 
available  at  10-30% 


Derivatives  X2  -  20 — *-30%  using  Cd  G-As^  for  4.5  -  6.2  pm 
Receiver 

Telescope:  1  meter  diameter,  0.2-12  pm  range 
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Table  5.3. 


Aerosol  Profiling  Experiment  Extracted  from  Reference  3d 


PROFILING  OF  TROPOSPHERIC  CLOUDS  AND  AEROSOLS 


Description:  Measure  the  presence,  geometrical  profiles  and  backseattering 
of  thin  clouds  and  aerosols  -  including  data  on  the  ir  optical  thickness.  Co¬ 
aligned  Earth  radiation  budget  sensor  would  measure  flux  and  albedo  change's; 
co-aligned  passive  temperature  and  humidity  sensors  would  allow  assess¬ 
ment  of  cloud  and  aerosol  effects  on  passive  profilers.  Data  would  docu¬ 
ment  the  locations  and  extent  of  aerosols  and  clouds  and  permit  correlation 
with  radiation  changes  detected  by  passive  sensors.  This  information  is 
essential  for  accurate  radiation  transfer  models.  The  existence  of  pollution 
and  dust  clouds  (e.g.  ,  Sahara  dust)  over  large  areas  could  be  monitored 
and  mapped. 

Implementation:  Standard  Nd-YAG  lidar  (IX,  2X)  with  1  mr  FOV.  2  PMT 
detectors  (20  percent  QE  at  530  nm;  2  percent  QK  at  1060  nnt).  Night  or 
day  operation.  Data  system  capable  of  generating  spatial  protile  ot  return 
with  height  resolution  of  about  150  m.  Some  simulations  are  available 
and  sample  results  are  shown  in  figure  (lidar  at  160  km,  0.1  .1 ,  Ravleian 
atmosphere  only)  for  expected  single -shot  errors.  Simulations  under 
development  to  include  tropospheric  aerosols,  higher  laser  energy  an  1 
multi-shot  averaging  show  greatly  improved  signal  to  noise.  System 
will  be  limited  by  signaL  levels  for  very  low  density  aerosols  and  clouds, 
and  signal  measurement  errors  will  be  limiting  factor.  Achieving  :  20 
percent  accuracy  in  conversion  to  optical  thickness  will  be  difticup  a  priori 
but  will  be  aided  by  2  backseattering  measurements.  Even  poorer  accuracy 
will  be  very  useful  in  most  applications. 

Feasibility:  1A  (530  nm),  IB  (1060  nm).  Experiment  design  information 

is  available,  and  a  moderate  amount  of  ground-based  field  data  is  available 
in  the  visible  region  of  the  spectrum.  The  requirement  for  good  range 
re  solution  at  very  low  signal  levels  increases  the  difficulty. 

heeded  Developments:  Optimization  of  the  1060  nm  svstem  i.-.  reurivd  atri 
additional  simulations  of  particular  ras<  s  are  nre  i,-:. 

Discussion:  At  present  no  other  svstem  exists  for  in«  asurme  ilv<  sna'ial 
distribution  of  subvisual  tropospheric  clouds  am:  a  r.vots.  -hier  ;:v.  nr. 
vital  tor  a  more  uuaut  iiativ,  understanding  >•  tin  .  auin’  ion  t  uu. 

’tic  operation  ot  mn  uitlv  ss-ui  passive  sct'svi.-s.  !r  addition.  moojtonnr 
u e  up. >•  sons  thi  u>  ;  'sols  and  'douus  nmo  :  s  m:.  •  i  .  >n  . 

!  rail.  T'D!  ‘  .t*  •  i  ‘ 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  review  set  out  to  answer  the  question  of  whether  there 
existed  in  the  published  literature  any  methods  for  remote  sensing  atmos¬ 
pheric  visibility.  Our  studies  have  led  us  to  the  conclusion  that  although 
there  exist  several  promising  techniques  of  remotely  monitoring  atmos¬ 
pheric  visibility  there  is  no  proven  technique  which  can  be  immediately 
put  to  use.  In  fact  only  two  or  three  methods  have  actually  been  subjected 
to  testing  using  real  satellite  observations.  Also  it  should  be  recalled  that 
these  methods  do  not  monitor  visibility  itself  but  rather  the  total  optical 
transmission  of  the  atmosphere.  A  separate  calculation  must  be  performed 
to  determine  visibility  for  the  actual  optical  system  and  environmental 
situation  under  consideration.  More  specifically,  the  three  concepts  which 
have  been  subjected  to  some  degree  of  experimental  verification  are: 

1.  Contrast  Reduction  {See  Section  4.  1). 

2.  Upwelling  spectral  radiance  (See  Sections  4.  2  and  4.  3). 

3.  Angle -integrated  Radiance  Difference  (See  Section  4.4). 

However,  more  experimental  measurements  will  be  needed  to  confirm  these 
tentative  conclusions  and  to  establish  the  sensitivity  and  range  of  applicability 
of  the  method.  An  important  adjunct  to  these  measurements  is  accurate 
measurement  of  the  ground-truth  values  used  for  comparison.  Because  of 
the  variability  of  atmospheric  conditions,  these  measurements  should  be 
carried  out  in  a  continuous  series  extending  for  a  minimum  of  2  to  3  months. 

In  the  course  of  this  review,  a  new  concept  for  measuring  the 
optical  depth  of  the  atmosphere  was  developed  that  makes  use  of  nighttime 
observation  of  natural  or  artificial  illumination  sources  on  the  earth  (See 
Section  4.  6).  Further  theoretical  and  numerical  calculations  should  be 
performed  to  test  the  possibilities  of  this  technique.  If  its  capabilities 
pass  these  additional  tests,  consideration  should  be  given  to  verifying  this 
method  using  existing  observing  satellites. 


Finally,  we  make  the  conjecture  that  any  operational  system 
would  depend  not  on  a  single  method  but  would  be  a  composite  of  several 
techniques  mutually  supporting  and  checking  the  others.  To  this  end,  we 
believe  that  as  many  as  possible  of  the  candidate  methods  should  be 
carried  forward  through  the  various  development  stages. 
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APPENDIX  A 


Derivation  of  Working  Equation  for  Simulated  Polar 
Nephelometer  Method  of  Section  4.4 

Figure  A1  shows  the  sun  at  zenith  angle  6  illuminating  the 

s 

atmosphere.  The  radiance  of  the  atmosphere  is  observed  at  two  directions 
symmetrically  located  about  the  zenith  at  the  same  angle.  The  differential 
equation  (DE)  governing  the  transport  of  radiation  along  the  observation 
path  considering  only  single  scattering  is 


-HU  +i3„)ds' 
Ja'  'aer  R 


-(0  +d„)B  +(d  cf,  ,  +  ,)H  e 

aer  R  li  aer  aer  1  R  R  1  o 


~f(6  +  0  )d  s 

o  aer  R 


+  BGie 


(A.  1) 


The  meaning  of  the  symbols  is  as  follows: 

Bli  ^  radiance  at  position  s  along  path  1  at  angle  0  to 
zenith,  w  cm  sr"1  1 

s  =  distance  along  path  1  measured  from  origin  at  surface,  km 

daer  r  scattering  coefficient  for  aerosol,  km"1 

=  Rayleigh  scattering  coefficient  for  air  molecules,  km  1 

^aer  1=  scattering  phase  function  for  aerosol  evaluated  for 
scattering  into  direction  1 

i  =  Rayleigh  scattering  phase  function  for  air  molecules 
evaluated  for  scattering  into  direction  1 

s'  -  distance  along  path  to  sun  measured  from  origin  at  surface,  km 

l 

Rpi  r  radiance  of  ground  surface  in  direction  of  path  1, 
w  cm"“  sr~ 

11  r  irradiance  of  sun  at  top  of  atmosphere,  w  cm 


The  differential  path  elements  ds  and  ds'  can  be  related  to  the  differential 
element  of  distance  in  the  zenith  direction  dx  as  follows 


ds  =  sec  Q.  dx  =  m.dx  ,  ds'  =  sec  6  dx  =  m  dx  (A.  2) 

11  s  s 

where  the  approximate  equality  between  the  secant  of  the  zenith  angle  and 
atmosphere  mass  m  in  the  direction  having  that  angle  has  been  used.  Assuming 
that  aerosol  absorption  is  small,  the  optical  depth  of  the  atmosphere  in  the 
vertical  direction  can  be  approximated  as 


/•* 

I  (3  +  3  )  dx  =  r  (x)  +  T  (x)  =  r  (x) 

/  aer  R  aer  R 

J  o 


(A.  3) 


f  +  0  )  dx  =  r 

/  aer  R  o 


+  T  =  T 

aer  oR  o 


(A. 4) 


for  the  entire  atmosphere.  Also  the  differential  optical  depth  is 


d  r  =  (0  +  3  )  dx 

aer  R 


Now  the  DE  for  path  1  can  be  rewritten  as 


-  rm  ( 3  3  )  B  +  m .  (  3  +S  tp.tH  e 

1  aer  +  R  li  i  aer  aerl  R  Rl  o 


-m  (  T  -  r(x)) 

s  o 


-m  r  (x) 

+  miBoi e  1 


(A.  1') 


A  similar  equation  for  path  2  follows  immediately. 


I 


Subtracting  the  two  equations  leads  to  a  DE  for  the  difference 
in  radiance  observed  in  the  two  directions  1  and  2. 


d  a  B 
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(A.  5) 


Now  the  last  term  on  the  r.h.  s.  is  made  to  vanish  by  invoking  the  assumption 
that  the  ground  radiance  is  symmetrical  about  the  zenith  angle.  The  factor  in 
parenthesis  in  the  second  term  on  the  r.h.  s.  is  rewritten  as  follows: 
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Thus,  the  DE  in  a  is 
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Then  after  dividing  through  by  {  8aer  +  3^)  and  transforming  to  the  optical 
depth  as  the  independent  variable,  the  DE  becomes 

3t 
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This  DE  can  only  be  integrated  analytically  if  the  term  in  the  square  brackets 


on  the  r.  h.  s.  is  independent  of  the  optical  depth.  This  is  equivalent  to 
assuming  both  that  the  ratio  +  anc*  P^ase  functions  are 

independent  of  optical  depth.  Under  these  assumptions,  the  DE  becomes 


d  msT 

Tr  “  B12  +  mi  1  B12  ”  mi  K  e 


(A.  9) 


The  first  order  linear  DE  equation  is  readily  solved  using 
Laplace  transforms  with  p  as  the  transform  variable  and  the  caret 
identifying  the  transformed  function. 


/s 

(p  +  rr^)  A  B 


12  p  -  m  +  A  ®12 
^  s 


(A.  10) 


Since  A  B, 


=  o  ,  the  equation  is 


mi  K 


12  (p  +  m- )  (p  -  m  ) 
1  « 


(A.  11) 


Applying  the  inverse  transform,  the  difference  in  radiance  is 

m.  m  r  -m.T 

A  B  „  (t)  =  - J -  K  (e  s  -  e  1  ) 

12  m.  +  m 

1  s 


(A.  12) 


-m  r  m  t  -m.T 


{  f  A*aerl2  U'*,R12"A4,aerl2>  1  6  }  (e 

-  °  '  aer  R 


m .  H  0 

_ I f\<b  + — 2 -  a4>  )1» 

m  +m  aer  12  0  +0  R12  aer  12 

i  s  aer  R 

-m  -  r  ) 
,  s  o 


-m  t  -m.T 

so  1  > 
e  ) 


When  evaluated  at  t  =  i.  e.  from  outside  the  atmosphere,  the  expression 


becomes 


m.H 
1  o 


12  o  m  +  m 
1  s 


fA4aerl2  +/3  +  73“  (A4>R12"  A*aerl2} 

aer  R 


-(m.  +  m  )  t 

it  1  s  o 

(  1  -  e  ) 


(A.  13) 


Now  make  the  following  algebraic  manipulations. 


A<t>  +  — — ■* —  (4*  -  A  <t>  )  = 

aerl2  aer  PR  R12  aerl2 


A  B  (  r  )  (m.  +  m  ) 
12  o  i  s 


H  m.  (1-e 
o  i 


-(mi  +  m>0 


(A.  14) 


(3  +  3  )  i<t  +  B  (A4>  -  A  4>  ) 

aer  R;  A  aerl2  pr'  R12  aerl2; 


A  B  (m.  +  m  )  (3  +  /3  ) 

12  i  s  aer  R 

-(m.  +  m  )  t 

H  m.(  1-e  1  S  °)  (A.  15) 

o  1 


If  the  phase  functions  4>aer  12  and  ^12  are  indePendent  of  altitude,  an 
integrationover  the  height  of  the  atmosphere  can  be  done  and  it  produces 


T  a  <t> 

o  aer  a*  r!2 


r  A  4> 
oR  R  12 


A  B  (m.  +m  )(t  +  r  ) 

12  1  so  aer  oR 


H  m.  (  1  -  e 
o  1 


-(m.  +  m  )  r 
l  s  o. 


(A.  16) 


The  Rayleigh  phase  function  is  always  independent  of  altitude  (neglecting 
depolarization).  For  the  aerosol  phase  function  to  be  independent  of  altitude 
requires  that  the  size  distribution  and  composition  be  independent  of  altitude 
which  in  general  is  not  true.  However,  in  the  boundary  layer  which  contains 
most  of  the  optical  thickness  of  the  atmosphere,  this  condition  of  altitude 
independence  is  approximately  true. 
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The  final  step  is  to  integrate  both  sides  of  Eq.  (A.  16)  over 
the  scattering  angle.  To  arrive  at  the  correct  expression  for  the  optical 
depth  of  the  atmosphere,  it  must  be  assumed  that  the  integral  over  the 
phase  function  in  the  back  scattering  region  can  be  neglected  compared 
to  the  forward  scattering  contribution.  The  consequence  is  that  the  method 
cannot  be  absolute  but  a  correction  or  calibration  factor  must  be  obtained 
from  a  separate  experiment. 

The  formal  solution  which  results  after  integrating  over  all 
scattering  angles  and  using  the  above  assumption  is 

7 r 

r  2  *  f  <t>  sin  0  d  6  =  r  (A.  17) 

o  aer  /  aer  1  o  *er 

•A) 


/TT 
O 


d  e  sin  e 


A  B.  „  (  t  +  T  )  (  m.  +  m  ) 

12  o  aer  o  R  i  s 

-(m.  +  m  )  r 

TT  /  1  1  SO 

H  m.  (  1  -  e  ) 

o  1 


T  A<t> 

R  R12 


which  is  the  working  equation  for  this  method.  The  analogous  working 
equation  for  the  ground  based  method  is 


T 

o  aer 


7r 

Jo 


d  Q  sin  6 


A  B.  „(  T  +  T  )  (m.  -  m  ) 

12  o  aer  o  R _ i _ s_ 

-  r  m  -  t  m 
TT  ,  os  o  1. 

H  m.(  e  -  e  ) 

o  1 


r  A 

R  R 12 


(A. 18) 
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APPENDIX  B 


Solution  of  Transcendental  Equation  for  ^  and  d 


In  Section  4.  5,  a  method  is  given  for  obtaining  the  vertical 
optical  thickness  from  measurements  of  the  angular  variation  of  the 
nighttime  intensity  of  a  city.  The  unattenuated  intensity  will  have  a  complex 
angular  variation  which  is  unknown.  The  analysis  of  Section  4.  5  assumed 
that  the  variation  had  only  a  constant  term  and  a  cosine  term  with  the 
relative  contribution  of  the  two  unknown.  A  single  transcendental  equation 
relating  the  vertical  optical  thickness  r  and  a  term  d  which  gives  the  pro¬ 
portion  of  the  cosine  variation  was  derived.  This  appendix  gives  a  numerical 
method  for  solving  this  equation  when  measurements  are  given. 

The  basic  equation  is 


j(Vw)  ,  (  !-<!<—  >1 


-  T  W 

o 


(B.  1) 


where  w  =  sec  0-1  and  w/  (w  +  1)  =  1  -  cos  0, 
An  easy  algebraic  manipulation  gives 


“  {j 1 


1  -  d  ( 


—  >J  } 


(B.  2) 


=  —  In  /  4-  [  1  -  d  ( 
o  w  l  ] 


w  +1 


(B.  3) 


It  is  also  easy  to  solve  for  d,  i.  e. 


J  W  +  1  r  *  .  TC  1  , 

d  *  — —  [ 1  -Je  1 


(B.4) 


Now  write  down  the  two  equations  for  r  using  two  measurements  j  and  j  . 

o  1  « 
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T 


(B.  5) 


—  In  i  i-  [  1  -  d  < - i-r  )  ]  l 

\  Uj  w1  +  1 

—  In  \  4-  [  1  -  d  ( - 2—  )  ]  l 

r2  *  32  W2  +  1  f 

Figure  B1  shows  that  when  plotted  as  a  function  of  d,  these  two  curves  inter¬ 
sect  and  give  a  common  solution  for  and  d. 

To  find  the  numberical  values  at  the  intersection,  rewrite  one  of 
the  equations  in  the  form  giving  d  explicitly  ,  say 


.  W2+1  r,  .  'ow2, 

d  “  “VT—  1  1  -  e  ’ 


(B.  7) 


Make  an  initial  estimate  of  where  the  solution  is  located.  Insert  the  value 

of  d  into  the  equation  for  r  .  Use  the  value  to  obtain  a  new  value  of  d. 

o 

Continuation  of  these  steps  produces  a  sequence  of  values  of  T  q 
and  d  which  converge  to  the  desired  solution.  A  geometric  diagram  of  the 
process  is  given  in  Figure  B2.  For  the  process  to  converge,  the  sequence 
must  be  to  find  from  function  1  and  d  from  function  2.  Convergence  is 
obtained  regardless  of  whether  the  initial  value  is  greater  than  or  less  than  the 
solution.  If  the  roles  of  1  and  2  were  reversed  in  this  example,  the  sequence 
of  values  would  diverge.  One  concludes  that  convergence  depends  on  the 
relative  slopes  of  the  two  curves.  Thus  if  divergence  occurs,  a  reversal  of 
the  equation  should  induce  convergence. 

This  process  is  readily  carried  out  with  a  small  programmable 
calculator.  All  of  the  results  given  in  Section  4.  5  were  obtained  with  a 
Hewlett-Packard  HP -55  calculator.  Table  B1  shows  the  rapid  convergent 
sequence  for  an  example  where  the  correct  solution  is  r  =  0.  1  and  d  *  0,  2. 
Values  of  and  j  correct  to  4  decimal  places  are  used.  Two  initial  values 
for  d  are  tried,  0  and  1,  which  are  the  extremes  of  that  variable. 
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Schematic  diagram  of  iterative  sohrion  of  liquations 
Ifi  and  1 17  . 


.1 


Table  B1  .  Examples  of  Converging  Sequence  Values  in  the 

Solution  of  the  Transcendental  Equation  for  and  d. 


CASE  A  CASE  B 


d 

To 

d 

To 

0.00000 

. 11054 

1.00000 

-. 11973 

. 18792 

. 10133 

.44951 

.06753 

. 19861 

. 10015 

.23770 

.09572 

. 19998 

. 10000 

. 20512 

.09943 

.20016 

.09998 

. 20082 

.09991 

. 20018 

.09998 

. 20026 

.09997 

.20019 

. 09998 

.20018 

. 09998 

APPENDIX  C 


First  Order  Downward  Scattering  of  Sunlight 


We  consider  sunlight  incident  on  a  clear  atmosphere  from  any 
zenith  angle  and  ask  the  question--for  those  photons  which  are  scattered 
in  the  atmosphere  before  reaching  the  ground,  what  is  the  probability  that 
the  first  scattering  event  will  direct  the  photon  in  a  downward  direction? 

Here  we  will  answer  the  question  at  wavelength  X  =  0.45^  simply 
because  we  have  the  necessary  tools  at  hand  to  do  this,  Deirmendjian's 
haze  model  L  (continental  haze)  is  used  to  represent  the  scattering  charac¬ 
teristics  of  the  aerosols  involved.  The  optical  depth  of  the  atmosphere 
when  the  zenith  angle  considered  is  zero  r  (o)  is  described  here  as  follows: 


r  (o) 


(h)  +  a  (h) 1  dh 
m 

(h)  dh 


(h)  dh  + 


(C.  1) 


where  a  (h)  is  the  aerosol  scattering  coeffieient  at  altitude  h 

cr^Oi)  is  the  molecular  scattering  coefficient  at  altitude  h 


The  contribution  of  the  molecular  scattering  may  be  written 


(h)  dh  «  8  cj 

m 


(o) 


T  (o) 
m 


(C.  2) 


-2 

and  at  X=  0. 45^  cr  =  2.72  x  10  or  t  (o)  =  0.217.  The  normalized 

m  m 

phase  function  for  molecular  scattering  of  unpolarized  light  <t>  (a)  is  given  by 


*(a)  = 


16  7T 


^  2  .  -1 
(1  +  cos  a)  ,  sr 


(C.  3) 


where  a  is  the  scattering  induced  angle  of  deviation  of  the  incident  ray. 


*  Height  of  homogeneous  atmosphere  =  8  km 
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This  implies  that  the  three  dimensional  uistribution  of  the  directional 
probability  of  scattering  has  a  center  of  symmetry.  Hence,  no  matter 
what  the  angle  of  incidence  on  the  atmosphere  of  light  the  probability 
that  light  scattered  by  molecules  will  be  deviated  into  a  downward  direction 
is  0.  5 

This  contribution  to  the  optical  depth  made  by  the  aerosols,  r  , 
depends  on  their  concentration.  If  the  optical  depth  r  is  known  or  assumed 
then 


r  (o)  =  t  (o)  -  0,  217  (C.4) 

a 

and,  more  generally,  if  0  is  the  zenith  angle  at  which  the  optical  depth  of 
the  atmosphere  r  (0)  is  evaluated  then 


r  (< p )  =  t  (o)  sec  0  =  t  (0)  -  0.  217  sec  <p  (C.  5) 

3  cl 


Assuming  no  absorption  in  the  atmosphere,  the  fraction  of  the 
light  incident  at  zenith  angle  0  which  is  scattered  before  reaching  the 
ground,  s,  is  given  by 


s  = 


(C.  6) 


The  fraction  scattered  by  molecules. 


s  is 
m 


s 

m 


*  (0) 

m 


(0) 


T  m  (°  )  0.  217 

r  ( O )  t  ( O ) 


The  fraction  scattered  by  aerosols,  s  ,  is 

’a 


r  (o) 

a _  _  t 

r  (o) 


(C.  7) 


(o)  -  0.  217 
r  (o) 


(C.  8) 


Normalized  phase  function,  <t>  (a) 
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Figure  C.  1.  Normalized  phase  function  for  unpolarized  light  at  0.45 
microns  as  given  by  Deirmendjian  (solid  curve)  and  calculated  by  fitted 
Equation  C.  9  (circled  points). 
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and  where  f  (a)  is  used  in  the  range  90°  <  a  <  180° 


90°  -  $ 

f<*>=  0  for  -  1 


and  f (o-) 


cos 


-1/90°  -  <p  \ 
' 180°  -  a  > 

180  ° 
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The  summing  approximation  illustrated  by  Equation  Cll  was  used  to  obtain 
values  of  R  (<p)  from  Equation  C12.  Values  calculated  are  presented  in 
column  1  of  Table  Cl  and  are  plotted  in  Figure  C2. 

The  second  column  of  Table  Cl  gives  the  fraction  of  molecular 
scattered  light  directed  in  a  downward  direction.  The  third  and  fourth 
columns  give  this  fraction  for  two  values  of  the  vertical  optical  depth  of  the 
atmosphere.  Equations  C7  and  C8  were  used  to  calculate  the  fraction  of  the 
total  scattering  due  to  aerosols  and  due  to  molecules. 
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Aerosol  Downward  Scattered  Fraction 


0  10 

20  30  40  50  60  70  80  90  \  ' 

Zenith  Angle  Degrees  J 

Figure  C. 2. 

First  order  downward  scattered  fraction  of  light  | 

scattered  at  0.45  microns  by  Model  L  aerosols. 

1 

' 

TABLE  Cl 


Zenith 

Angle 

Aerosol 

Downward 

Fraction 

Molecular 

Downward 

Fraction 

R  (o) 

t  =  0.  6 
Downward 
Fraction 

r  =  0.4 

Downward 

Fraction 

0 

.  947 

.  50 

.785 

.705 

10 

.  947 

.  50 

.785 

.705 

20 

.  946 

.  50 

.785 

.704 

30 

.  943 

.50 

.  783 

.703 

40 

.  937 

.  50 

.779 

.700 

50 

.  924 

.  50 

.  771 

.  694 

60 

.  898 

.50 

.754 

.  682 

66 

.  871 

.  50 

.737 

.  640 

70 

.  844 

.50 

.720 

.  657 

76 

.785 

.  50 

.  682 

.  630 

80 

.729 

.  50 

.  646 

.  605 

82 

.  693 

.  50 

.  623 

.  588 

84 

.  658 

.  50 

.  601 

.  572 

86 

.  607 

.  50 

.  568 

.  549 

88 

.  556 

.  50 

.  536 

.  526 

90 

.  500 

.  50 

.  500 

.  500 

First  Order  Downward  Scattered  Fraction  of  Scattered  Light, 
R  (  <M  at  \  ~  0. 45p  for  Haze  Model  L,  for  Molecular 
Scattering  and  for  Two  Values  of  Vertical  Optical  Depth  of  the 
Atmosphere. 
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